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Abstract 

The influence of hole-hole propagation in addition to the conventional particle- 
particle propagation, on the energy per nucleon and the momentum distribution is 
investigated. The results are compared to the Brueckner-Hartree-Fock (BHF) cal- 
culations with a continuous choice and conventional choice for the single-particle 
spectrum. The Bethe-Goldstone equation has been solved using realistic NN inter- 
actions. Also, the structure of nucleon self-energy in nuclear matter is evaluated. All 
the self-energies are calculated self-consistently. Starting from the BHF approxima- 
tion without the usual angle-average approximation, the effects of hole-hole contri- 
butions and a self-consistent treatment within the framework of the Green function 
approach are investigated. Using the self-consistent self-energy, the hole and particle 
self-consistent spectral functions including the particle-particle and hole-hole ladder 
contributions in nuclear matter are calculated using realistic NN interactions. We 
found that, the difference in binding energy between both results, i.e. BHF and 
self-consistent Green function, is not large. This explains why is the BHF ignored 
the 2hlp contribution. 
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1 Introduction 



It is one of the fundamental issue in nuclear physics to evaluate the nuclear 
matter binding energy and saturation properties, starting from a realistic 
nucleon-nucleon (NN) interaction without any adjustment of a free parameter 
[1,2,3]. The solution of this problem is of great interest, because the under- 
standing of the binding energy of nuclei is one of the basic problems of nuclear 
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physics and physics in general. Because of the presence of strongly repulsive 
components in the short-range part of the NN interaction [4,5,6], the NN 
scattering correlation has a predominant importance in the nuclear matter 
calculation. The solution of this problem is also very important because the 
nuclear many-body problem is one of the most challenging testing grounds 
for many-body theories of quantum systems. Last not least, however, a micro- 
scopic understanding of the binding energy of normal nuclei, which is based 
on a realistic model for the NN interaction, shall also provide a reliable pre- 
diction for the equation of state of nuclear matter at densities beyond the 
saturation density of normal nuclei, which is required for the study of astro- 
physical objects like the explosion of supernova and the structure of neutron 
stars. 

Recently, Various approximation schemes have been developed to describe the 
correlations which are induced into the many-nucleon wave function by the 
strong short-range and tensor components of such a realistic NN force. For a 
recent review on such methods see, e.g., Refs. [7,8]. The resulting correlation 
functions are the key input for the study of the absorption of real and virtual 
photons by pairs of nucleons as it is observed in (7, NN) or (e, e' NN) 
reactions. One of the most popular approximation schemes, which is frequently 
used in nuclear physics, is the Brueckner-Bethe-Goldstone (BBG) theory [1] 
which has been one of the tools for solving the many-body problem of nuclear 
matter for many years [8,9,10,11]. The results of BHF calculations depend on 
the choice of the single particle potential, U(k). In the conventional choice, 
U = for k > hp-, and U is the self-consistent BHF potential for k < kp. The 
alternative "continuous" choice has been proposed in [12] for which U is again 
the self-consistent BHF potential, but it extends to k > kp. This continuous 
choice leads to an enhancement of correlation effects in the medium and tends 
to predict larger binding energies for nuclear matter than the conventional 
choice. 

In this paper, special attention will be paid to study the self-consistent self- 
energy. Note, that all the self-energies in this paper will be calculated self- 
consistently. The self-energy of a nucleon contains all the information neces- 
sary to obtain occupation probabilities, quasi-particle strength and broaden- 
ing features which alternatively can be visualized in terms of hole and particle 
spectral functions. 

It is the purpose of this paper to investigate whether an improven and more 
consistent treatment of the many-body system, still based on nonrelativistic 
dynamics and two-body forces, can provide new information on the saturation 
problem of nuclear matter. A self-consistent treatment of the Green's functions 
formalism leads to the so called self-consistent Green's function (SCGF) the- 
ory [13,14,15,16,17,18,19]. Attempts have been made to employ the technique 
of a self-consistent evaluation of Green's functions [20,21] to the solution of 
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the nuclear many-body problem. This method offers various advantages: (i) 
The single-particle Green's function contains detailed information about the 
spectral function, i.e. the distribution of single-particle strength, to be ob- 
served in nucleon knock-out experiments, as a function of missing energy and 
momentum, (ii) The method can be extended to finite temperatures [5,22], 
a feature which is of interest for the study of the nuclear properties in as- 
trophysical environments, (iii) The BHF approximation, the approximation 
to the hole-line expansion which is commonly used, can be considered as a 
specific approximation within this scheme. Attempts have been made to start 
from the BHF approximation and include the effects of the hole-hole scatter- 
ing terms in a perturbative way [23,24]. There is a good reason to believe that 
a more complete treatment of the nucleon self-energy in the full framework 
of Green-function theory will ultimately help resolve long standing questions 
about nuclear saturation. Via this study we will be able to give an answer to a 
fundamental question regarding reliability of the BHF theory. Application of 
SCGF theory has proven useful for understanding the location and quantity 
of single-particle strength found in the electron-scattering experiments. The 
focus of this paper is to apply SCGF theory to nuclear matter to study the 
nucleon properties. 

After this introduction we will discuss some features of the self-consistent BHF 
approach for the self-consistent self-energy in Sec. 2. The effect of the hole- hole 
terms in the self-consistent self-energy and a self-consistent Green-function will 
be presented in Sec. 3. The self-consistent spectral function will be presented 
in Sec. 4. The main conclusions are summarized in the final section. 



2 BHF approximation 



One of the central equations to be solved in the BHF approximation is the 
Bethe-Goldstone equation. It is given by 



< k' x k' 2 \G(uj)\k 1 k 2 > = < k[k 2 \v\kik2 > + ^2 < k^k^v^k^ > 

x Qi i 3 ; u h l, < hhio^ihh > , (i) 

u - E(k 3 , fc 4 ) 

where v is the bare potential, u denotes the starting energy, Q is the Pauli 
operator which requires the nucleon momenta to be outside the Fermi sea, E 
is the sum of the two single-particle energies inside nuclear matter given by 

E(ki,k 2 )=e kl + e ka , (2) 
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with the single particle energies 

e k = -^ + ReX BHF (k,u; = e k ), (3) 
2m 

where the self-energy of a nucleon in nuclear matter with momentum k is 
determined by the self-consistent equation 

T, BHF (k,uj) = f d 3 k' < kk'\G{e k + e w )\kk' > n (k') , (4) 



with the occupation probability of a free Fermi gas with a Fermi momentum 
k F 

{ 1 for \k'\ < k F 
n (k') = { 'J" . (5) 

[0 for |Jfe'| > kp 

In the BHF approximation, (1), (3), and (4) are solved self-consistently. 

The definition of the single-particle potential or the self-energy can be ex- 
tended beyond the BHF approximation. In order to discuss such an extension 
we want to consider various contributions, which are represented by the dia- 
grams or series of diagrams displayed in Fig. 1. These terms are summed up 
in a self-energy as can be written schematically by 



E(k,u) = j: HF (k) + Aj: 2plh (k , uj) + AT l 2hlp (k, uj) 

= E BHF (k J uj) + AE 2hlp (k J uj) 1 (6) 

where the contribution due to two-particle one-hole terms AJ2 2plh is complex 
for uj > Bp while the two-hole one-particle term A£ 2/!lp is complex for u < ep. 
The imaginary part of the BHF self-energy vanishes for uj < ep. 

Real and imaginary parts of the self-energy are related to each other by a 
dispersion relation of the form [25,26] 

ReE BHF (k t uj) = Z HF (k) + -J Im Z BHF (k,uj') dujl 

71 J UJ 1 — UJ 



— oo 

oo 



A 261n/1 . 1 f ImAZ 2hl P(k,uj') , . 
ReAE 2hlp (k, uj ) = - / y -^—^duj', 8 

71 J Uj' — UJ 

— oo 

The self-energy of the Hartree-Fock term is real and independent of uj. It is 
given as [27] 

T J HF {k) = < kk' | V | kk' > n(k'), (9) 

k> 

where n(k') is the single particle occupation probability. The importance of 
the Hartree-Fock term is that it represents the simplest solution to the many 
body problem. 
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The energy dependence of the real part of the BHF self-consistent self-energy 
(see (4)) is visualized in Figs. 2 and 3, for nuclear matter with Fermi momen- 
tum kp = 1.36 fm _1 . It is evaluated for the CD-Bonn [28] and Bonn C potential 
[29], for various momenta k. This real part displays a pronounced minimum 
at energies around the Fermi energy. These results have been obtained after 
establishing self-consistency for the single-particle spectrum according to (3) 

Figs. 4 and 5 show the imaginary part of the BHF self-consistent self-energy 
as a function of u evaluated for the CD-Bonn (Fig. 4) and Bonn C (Fig. 5) 
potentials, respectively, for different values of momenta k. We see that the 
imaginary part is identical zero for energies u less than — Sp, as can be 
seen from (3), and yield non-negligible values up to very high energies. Al- 
though not explicitly shown in Figs. 4 and 5, lm'E BHF (k, u) vanishes at high 
energies which is a consequence of the vanishing of the imaginary G matrix 
elements at very large values of the energy parameter. This result is general 
for any soft-core interaction, although the softness of the core directly influ- 
ences the energy for which the matrix elements become negligible. Considering 
the results for the imaginary part of the self-consistent self-energy, which are 
displayed in Figs. 4 and 5, it is clear from the dispersion relation (7) that the 
real part of the self-energy is identical to the corresponding HF single-particle 
potential H HF (k) in the limit uj — > —oo. This real part gets more attractive 
with increasing u until one reaches values of u at which the imaginary part 
is different from zero. The self-energy turns less attractive at higher energies, 
which leads to a pronounced minimum at energies u slightly above the Fermi 
energy. 

The conventional choice has been to ignore self-energy contributions for the 
particle states completely and approximate the energies by the kinetic energy 
only. This conventional choice for the single-particle spectrum, however, is 
not very appealing as it leads to a gab at the Fermi surface: the propagator 
for single-particle states with momenta below the Fermi momentum kp is 
described in terms of a bound single-particle energy while the corresponding 
spectrum for the particle states starts at the kinetic energy for the momentum 
kp, so the potential in this case is given by [30] 



Alternatively, a continuous choice for the potential can be made [29] by defin- 
ing 



k'<k F 

for all states k below and above the Fermi surface. This leads to a spectrum 
that is continuous at the Fermi momentum. Where | kk' >a = \ kk' >-| k'k >. 




(10) 



U(k) = Re < kk ' I G ( u = e k + e*0 I kk' > A , 



(11) 
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The self-consistent single-particle potential can be written as follows 



BHF 



U B HF(k) = Re J2(k,u = e k ) . 



(12) 



Fig. 6 displays the single-particle self-consistent potentials as a function of mo- 
mentum at different densities using the CD-Bonn potential. The single-particle 
potential at kp=1.36 fm ~ 1 shows a significant deviation from a parabolic 
shape in particular at momenta slightly above the Fermi momentum, see also 
[31,32,33]. It is obvious that such a deviation tends to provide more attractive 
matrix elements of G in evaluating the self-energy for hole states according to 
(4), which leads to more binding energy. At high Fermi momentum (kp = 2.04 
fm -1 and kp = 2.72 fm -1 ) this deviation is disappear. 



3 Self-consistent Green functions 

In the self-consistent Green's function (SCGF) formalisms, the self-energy 
is the key quantity to determine the one-body Green's function. One can 
calculate the contribution of the hole-hole terms to the self-energy in a kind 
of perturbative way [34] 



We assume a single-particle spectrum e k which is identical to the self-consistent 
BHF spectrum, but shifted by a constant C±, which ensures the self-consistency 
for k = kp 



= -^ + Re [E BHF (k Fl u = e kp ) + AZ 2hl ?(k Fl u = e kw )\ . (14) 



This shifted single-particle spectrum is also used in the Bethe-Goldstone equa- 
tion. 

Results for the two-hole one-particle contribution to the self-consistent self- 
energy, AS 2Wp , are displayed in Figs. 7 and 8, for various momenta k. Fig. 8 
illustrates the self-consistent solutions of the imaginary part of the self-energy, 
including only 2hlp. The imaginary part of A£ 2Mp is different from zero only 
for energies uo below the Fermi energy. This implies that Im Tt 2hlp (k,u) = 0, 
for uj = Bp. From (8) one also observes, noting that the imaginary part of 
AE 2Wp (fc, u) is positive, that for k = kp the on-shell value of the real part 




(13) 




BHF 



+ G 



i 
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(u = ef) must be positive. It turns out that this holds for all single particle 
momenta below hp- The conservation of the total momentum in the two- 
nucleon of the G-matrix in (13), hi + /12 = k + p, leads to a minimal value 
of uj at which this imaginary part is different from zero. Due to these lim- 
itations the imaginary part integrated over all energies is much smaller for 
AS 2/llp than for T, BHF . The real part of A£ 2Wp is related to the imaginary 
part by a dispersion relation similar to the one of (8), connecting the imagi- 
nary part of Y] BHF with the particle-particle ladder contributions to the real 
part of Tj BHF . Since the imaginary part of AT, 2hlp is significantly smaller than 
the one of I] BHF , the same is true also for the corresponding real part. The 
contribution of AS 2Wp (fc, w) term is considerably smaller than the contribu- 
tion of Y* BHF (k,u). At first sight this seems to confirm the notation that the 
effects of hh-propagation are considerably smaller than those corresponding 
to pp-propagation present in G-matrices. 

The Lehmann representation of the Green's function for infinite nuclear matter 
is given by [35] 

sM = + J^W^L, (15) 

J LO — UJ — IT) J UJ — UJ + IT] 

— oo sf 

where Sh( P )(k,uj) is the hole (particle) spectral function and we will discuss it 
in detail in the next section. 

The formal solution of Dyson's equation is particularly simple in an infinite 

system, 



1 

u-£-UkM 



g{k,uj) = p — -. (16) 



In this approximation the solution of the Dyson equation leads only to a shift 
in the single particle energy according to (3). The Green's function expression 
for the binding energy is given by [36,31] 

E _ f<Pkfl* OB du>S h (k,u)±(£ : + u) 

A Jd 3 k n{k) " 1 ' 



Results for the binding energy per nucleon are displayed in Table 1. In this ta- 
ble we have compared results for the binding energy of nuclear matter using the 
CD-Bonn potential from different prescriptions; (i) Hartree-Fock, (ii) BHF (us- 
ing continuous choice), and (iii) BHF+2hlp (using conventional choice where 
the self-energy is calculated self consistently). The continuous choice leads 
to an enhancement of correlation effects in the medium and tends to predict 
larger binding energies for nuclear matter than the conventional choice. The 
conventional choice has been to ignore self-energy contribution for the particle 
states completely and approximate the energies by the kinetic energy only, see 
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(10). In nuclear matter one observes, that the inclusion of the 2hlp terms in 
the self-energy leads to less attractive quasiparticle energies, but more bind- 
ing energy. The gain in binding energy due to the 2hlp components in the 
self-energy is about 0.5 MeV per nucleon in nuclear matter at saturation den- 
sity [31]. We found that, no large difference in binding energy between both 
results, i.e. BHF and BHF+2hlp. 



4 Self-consistent spectral function 



The physical meaning of the spectral functions is rather simple. The hole spec- 
tral function Sh{k, uj) is the probability of removing a particle with momentum 
k from the target system of A particles leaving the resulting (A — 1) system 
with an energy E A ~ X = E A — uj, where E A is the ground state energy of the 
target. Analogously, the particle spectral function Sp(k,uj) is the probability 
of adding a particle with momentum k and leaving the resulting (A + l) system 
with an excitation energy uj measured with respect to the ground state of the 
A system, i.e., uj = E A+1 — E A . The spectral functions are obtained from the 
Dyson equation [25]. For the hole spectral function this yields the relation to 
the self-energy given by 



Sh(k,u) = —Ira g(k,uj) 

71 



Im Tj(k, uj) 



tt (uj - k 2 /2m- ReE(k,uj)) 2 + (/m E(fc,w)) 2 ' v ' 

The importance of the particle spectral function is to exhibit where the missing 
single-particle strength is located in energy. It is obtained by solving the Dyson 
equation and is related to the self-energy by 

a (k \ _ _1 ImZ(k,u;) 

p{ ,LJ> 7r(uj-k 2 /2m-ReZ(k,uj)) 2 + (Im?3(k,uj)) 2 ' 1 ] 



In Figs. 9 and 10 a typical example of a hole spectral functions, for three 
different momenta in nuclear matter at kp = 1.36 fm" 1 as a function of energy. 
The self-consistent spectral function (SCSF) calculation is performed for the 
CD-Bonn and Bonn C potentials. The self-consistent particle spectral function 
displays a peak at e qp because of the vanishing term in the denominator of 
(18). Notice that as k — > kp this peak becomes extremely sharp due to the 
vanishing Im T,(k,uj) in (18). 

In Figs. 11 and 12 a typical example of a particle spectral functions, for three 
different momenta in nuclear matter at kp = 1.36 fm _1 as a function of energy. 
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The SCSF calculation is performed for the CD-Bonn and Bonn C potentials. 
For momenta larger than kp a quasi-particle peak, which, as usual, broadens 
with increasing momentum, can be observed on top of the same high energy 
tail. The results therefore display a common, essentially momentum indepen- 
dent, high-energy tail. The location of single particle strength at high energy 
simply means that the interaction has sufficiently large matrix elements. 

The quasi-particle contribution can be isolated from the exact single particle 
propagator for k close to kp by expanding the self-energy around the Fermi en- 
ergy and employing the fact that the imaginary part of the self-energy behaves 
like (u> — Ef) 2 [37]. The quasi-particle energy is given by 

4 P co F (k) = ^ + ReE{k,ef caF (k)), (20) 

which coincides with (3). With these ingredients one obtains in the limit 
k — > kp a 5-function contribution to the single particle propagator. Since 
the imaginary part is typically small around sp, this quasi-particle behavior 
is also obtained for other momenta close to kp. 

Fig. 13 displays the single-particle energy e{k) (see (10) and (11)) as a function 
of the momentum, at kp— 1.2 fm -1 using CD-Bonn potential. In this figure 
we displayed the single particle energies for continuous choice, single particle 
energies below kp for SCGF using the conventional choice and the quasi par- 
ticle energy (20). To obtain correct sum rules for momenta close to the Fermi 
momentum kp, the quasi-particle approximation is used to isolate this part of 
the propagator explicity. The quasi-particle strength z(k) is given as [12] 

z(k) = \i-( dRe ^ k > u) ) r\ (2i) 

where the quasi-particle strength in nuclear matter at kp = 1.36 fm -1 using 
various models for the NN interaction are shown in Fig. 14. For comparison, 
results from the full many-body calculation (at kp = 1.33 fm _1 ) of Benhar et 
al. [38] (dots) are also shown in Fig. 14. 

Integrating the hole spectral function Sh{k,u) over all the accessible excited 
states of the (A-l) system, one obtains the occupation probability of the sin- 
gle particle momentum k in the correlated ground state. The single nucleon 
momentum distribution in nuclear matter is defined as 

n(k) — J Sh(k, u)du) — 1 — J Sp(k,u)duj. (22) 

— oo oo 

The momentum distribution is calculated self-consistent. Results are displayed 
in Fig. 15 for the CD-Bonn and Bonn C potentials at kp = 1.36 fm -1 . The 
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experimental data are taken from [39]. The self-consistent occupation proba- 
bilities for nuclear matter are in good agreement with experimental data. 



5 Conclusion 

In this paper the nucleon properties in the nuclear medium have been studied 
using a self-consistent self-energy. The effects of hole-hole contributions and a 
self-consistent treatment within the framework of the Green function approach 
are investigated. It is observed that the effects of hole-hole contributions in 
binding energies are not large, this represent a success for the BHF theory. 

Miither and his collaborators [31] and Baldo and Fiasconaro [33] shown, that 
the parabolic approximation for the single particle potential U(k) in the self- 
consistent Brueckner scheme introduces an uncertainty of 1-2 MeV near the 
saturation density, and therefore it cannot be used in accurate calculations. 
The full momentum dependence has to be retained, which prevents the use 
of a constant effective mass approximation. In this paper, we have explained 
that at high densities the parabolic approximation is valid. 

The spectral functions are calculated from the momentum and energy de- 
pendent self-consistent self-energy by solving the Dyson equation. We found 
that the nuclear matter self-consistent occupation probabilities are in good 
agreement with the experimental measurements. 
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Figure 1. Diagrams contributing to the nucleon self-energy. The wavy lines indicate 
the Brueckner G-matrix. (a) The left diagram is the standard Brueckner approxima- 
tion for the nucleon self-energy (4). (b) This diagram takes into account transitions 
to one-particle two-hole states. It contributes to both real and imaginary parts of 
the self-energy in the considered energy and momentum range. 




Energy a [MeV] 

Figure 2. The real part of the BHF self-consistent self-energy as a function of the 
energy uj (see (4)) for symmetric nuclear matter with Fermi momentum kp = 1.36 
fm _1 evaluated for the CD-Bonn interaction, calculated for various momenta k. 
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Figure 3. The real part of the BHF self-consistent self-energy as a function of the 
co (see (4)) for nuclear matter with Fermi momentum kp = 1.36 fm -1 evaluated for 
the Bonn C potential, for various momenta k. 




-500 500 1000 

0) [MeV] 

Figure 4. The imaginary part of the BHF self-consistent self-energy as a function of 
co (see (4)) for nuclear matter with Fermi momentum kp = 1.36 fm _1 calculated 
for the CD-Bonn interaction, calculated for various momenta k. 
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-500 500 1000 

0) [MeV] 

Figure 5. The imaginary part of the BHF self-consistent self-energy for nuclear mat- 
ter with Fermi momentum Hf = 1.36 fm _1 calculated for the Bonn C interaction, 
calculated for various momenta k. 
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Figure 6. The single-particle self-consistent potentials Ubhf (see (12)) as a function 
of momentum. Results are displayed at kp = 1.36 fm -1 , kp = 2.04 fm" 1 and 
kp = 2.72 fm" 1 using the CD-Bonn potential. 
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Figure 7. The real part of the 2hlp contribution to the self-consistent self-energy as 
a function of u (see (13)) evaluated for the CD-Bonn potential assuming kp = 1-36 
fm _1 for various momenta k. 




Figure 8. The imaginary part of the 2hlp contribution to the self-consistent 
self-energy as a function of lo (see (13)) evaluated for the CD-Bonn potential as- 
suming kp = 1.36 fm _1 for various momenta k. 
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Figure 9. Hole self-consistent spectral functions Sh{k, u>) as a function of u for three 
different momenta in nuclear matter at Icf = 1-36 fm" 1 assuming the CD-Bonn 
potential. Note the narrowing of the peak in the spectral function when k gets 
closer to hp. 
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Figure 10. Hole self-consistent spectral functions for three different momenta below 
kF in nuclear matter at kp = 1.36 fm" 1 using the Bonn C potential. 
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Figure 11. Particle self-consistent spectral functions S p (k,uj) as a function of en- 
ergy for various momenta showing a broadening quasi particle peak with increasing 
momentum. The data have been obtained for nuclear matter using the CD-Bonn 
potential with a Fermi momentum kp = 1-36 fm _1 and using gab choice. 
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Figure 12. Particle self-consistent spectral functions for k = 1.25, 1.5, and 2.0 kp. 
The results have been obtained for nuclear matter with a Fermi momentum kp = 
1.36 fm" 1 and using gab choice, assuming the Bonn C potential. 
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Figure 13. The single-particle energy e(k) (see (10) and (11)) as a function of the 
momentum at fcp=1.2 fm" 1 using the CD-Bonn potential. For the single particle 
energy for BHF using continuous choice, single particle energy below /cp for SCGF 
and quasi particle energy (20). 
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Figure 14. Quasi particle strength in nuclear matter at /cp=1.36 fm _1 (see (21)) 
using various models for the NN interaction. For comparison, results from the full 
many-body calculation of Benhar et al. [38] (dots) at kp= 1.33 fm _1 are also shown. 
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Figure 15. Momentum distribution in nuclear matter calculated from the spectral 
function. Results are given for the CD-Bonn and Bonn C potential at k F = 1.36 
fm . The experimental data are from [39]. 



Table 1 

Energy per nucleon for nuclear matter considering three different densities, kp = 
1.20, 1.36 and 1.60 fm , assuming the CD-Bonn potential. Results are displayed 
for HF, BHF with a continuous choice and using the exact Pauli operator [31], 
and the BHF+2hlp using conventional choice where the self-energy is calculated 
self-consistently. All energies are given in MeV per nucleon. 



k F [fm - 1 ] 


HF 


BHF, exact [31] 


(BHF+2hlp) 


1.20 


2.88 


-15.39 


-15.34 


CD-Bonn 1.36 


4.61 


-18.83 


-18.99 


1.60 


10.31 


-22.86 


-22.78 
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